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ABSTRACT 



Outflows generated by protostars heavily affect the kinematics and chemistry of the hosting molecular cloud through strong shocks 
that enhance the abundance of some molecules. LI 157 is the prototype of chemically active outflows, and a strong shock, called Bl, 
is taking place in its blue lobe between the precessing jet and the hosting cloud. We present the Herschel-PACS 55-210 /jm spectra of 
the L1157-B1 shock, showing emission lines from CO, H2O, OH, and [Oi]. The spatial resolution of the PACS spectrometer allows 
us to map the warm gas traced by far-infrared (FIR) lines with unprecedented detail. The rotational diagram of the high-y^p CO lines 
indicates high-excitation conditions (T^^ ^ 210 ± 10 K). We used a radiative transfer code to model the hot CO gas emission observed 
with PACS and in the CO (13-12) and (10-9) lines measured by Herschel-HIFI. We derive 200< ry„ <800 K and n > 10^ cm"'. The 
CO emission comes from a region of about 7" located at the rear of the bow shock where the [O i] and OH emission also originate. 
Comparison with shock models shows that the bright [O i] and OH emissions trace a dissociative J-type shock, which is also supported 
by a previous detection of [Fell] at the same position. The inferred mass-flux is consistent with the "reverse" shock where the jet is 
impacting on the LI 157-Bl bow shock. The same shock may contribute significantly to the high-/„p CO emission. 
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1. Introduction 

Shocks from protostellar outflows heat and compress the sur- 
rounding medium, triggering chemical processes that enrich the 
che mical composit i on of the environment of young stars (see 
e.g. lBachilleretal. I I2OOI . Although outflows and shocks have 
been studied since the eighties, we are still far from a clear un- 
derstanding of the interaction between the jet/outflow and the 
ambient cloud. In particular, studies based on the Infrared Space 
Observatory (ISO) data have shown that the atomic and molec- 
ular lines emitted in the FIR range significantly contribute to the 
cooling of the shocked material in outflows. However, because 
of the large ISO beam (80") only information averaged along a 
large part of the outflow lobe could be inferred, and a mixture of 
both continuous C-type and dissociative J-type shocks has been 
suggested to be responsible for the excitation of the observ ed 
FIR lines (e.g. lGiannini et al" l l200lHFlower & Gusdorfll2009b . 

In the outflow powered b y the L1157-mm Cl ass proto- 
star at a distance of 250 pc (iLoonev et al.L l2007l) . more than 
20 species have been detected in the sout hern, blue-shifted lobe 
jBachiUer et al.L 120011; lArce et al.L l2008h . The brightest shock, 
Bl, has been extensively observed in several molecular species 
at miflimetre (Gueth et al., 1996, 1998; Benedettini et al., 2002 



ICodella et all l2009h. near-infrared dNisini et al.L 
FIR wavelengths (' Giannini et alP 200r; Ni sini et al 



2010ah and 



2010b ') 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



re- 
vealing a complex structure with various shock tracers peak- 
ing at different positions. Previous modelling concluded that a 
(possibly non-stationar y) C-type shock is taking place in Bl 
(Gusdorf et al., 2008; N eufeld et a l., 2009; Viti et al., 2011). 

Because of its chemical rich ness, LI 157-Bl was se lected 
to be observed with the H IFI (Ide Graauw et al.L I201Q) and 
PACS_JPpglitsch e t all 1201 Ol) spectrometers on board Herschel 
('Pilbrattet al., 20 1(?), as part of the Chemical Herschel Su rvey of 
Star forming regions (CHESS) Key Program (Ceccarelli et all 
2010). The first results obtained by HIFI confirm the chemical 
richness of Bl (ICodella et al .l ^20 10*), including the first detec- 
tion of HCl in shocked gas (Codella e t al.Ll201ll) . They also re- 
veal a high-velocity gas with higher H2O abundance and excita- 
tion conditions than the bu lk of the emiss ion, and tracing a more 
compact region within Bl (iLefloch et al.Ll2010h . 

In this letter we present the PACS spectra of LI 157-Bl, fo- 
cusing our analysis on the high-Zup CO, OH and [O i] lines with 
the aim of investigating the high-excitation gas component. The 
H2O lines will be analysed in a forthcoming paper. 



2. Observations and data reduction 

The range spectroscopy mode of the PACS instrument was used 
to obtain the full range spectrum towards LI 157-Bl (RA(J2000) 
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Table 1. List of CO, [O i] and OH lines detected in spatial pix- 
els centred at offset (-5",7") and (0",0")- For lines covered in 
the two PACS scans, both measures are reported. Quoted errors 
(+lcr) include only statistical noise. Upper limits are 3cr. 
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"This line is affected by a large flux calibration error and it is not 
considered in the analysis. 



= 20''39™10;2, Dec(J2000) = +68°01'ia.'5). Two observations 
were carried out on May 25, 2010, covering the spectral ranges 
55-95. 2//m and 101.2-210 jim. The resolving power ranges be- 
tween 1000 and 4000, depending on wavelength and order The 
PACS field of view (FOV) of 47"x 47" is composed of 5x5 spa- 
tial pixels (spaxels) providing a spatial sampling of 9'.'4/spaxel, 
with a position angle of 61'?83. Our observations were carried 
out in staring mode, i.e. our maps are not Nyquist-sampled with 
respect to the point spread function (PSF). The data were re- 
duced with the HIPI4j v5.0 package. The flux calibration was 
performed by calculating the ratio between the source flux and 
the telescope background, and subsequently multiplying by the 
spectrum of Neptune. This procedure gives a flux accuracy of 
about 10% shortwards of 190 jum, while at longer wavelengths 
the flux calibration is much more uncertain, which caused us 
to exclude the CO (13-12) line at 200.25 fim in the analysis. 
Post-pipeline reduction steps were performed to extract the spec- 
trum of each of the 25 spaxels, to calculate the line fluxes with a 
Gaussian fitting of the line profile and to produce line intensity 
maps. 

Table [T] lists the line fluxes of the analysed species, namely 
CO, OH and [O i], detected at the centre of the FOV (the nominal 
position of Bl) and towards the PACS emission peak, located 
at offset (-5", 7") with respect to the centre of the FOV. The 
spectra of these two spaxels are shown in Figs. lA.l] and lA.2l The 
spectral range between 161.5 fim and 190.2 jum was observed 
in the two observations; the fluxes of the lines detected in both 
scans generally agree within the errors and in any case never 
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Fig. 1. Maps of the CO fines from (14-13) to (22-21). First con- 
tour and contour steps are 9x10""* W m^^/spaxel, corresponding 
to 3 cr. Crosses mark the centre of the spaxels of the PACS FOV. 
Circles represent the half -power beam width. 



differ by more than 20%, therefore we consider the averaged 
value in the following analysis. 



3. Results 

3.1. Morphology 

Figure [1] presents the maps of the CO lines from (14-13) to 
(22-21). The emission peaks towards spaxel (-5", 7"), with an 
elongation towards south-east. The deconvolved full with at half 
maximum (FWHM) size of the CO peak, estimated by sub- 
tracting in quadrature the FWHM of the PSF, ranges between 
5"and 10" depending on wavelength, with a mean value of 7", 
hence marginally resolved in these observations. The emission 
from [O i] and OH peaks in the same spaxel as high-Zup CO. In 
Fig.|2]we show the CO (14-13) and [Oi] 63 yum PACS maps 
superimpose d on tracers of postshock gas such as Si O (2-1 
(iGuefli et al.L [19981) and CH3CN (Sk-Ik) dCodella et al.L I2OO9I 
observed with the Plateau de Bure at a resolution of ~ 3". While 
both SiO and CH3CN trace an extende d bow shock at the tip of 
the outflow cavity traced by CO (1-0) dGueth et al .11 1996 1). it is 
clear that high-Zup CO, OH and [O i] are peaking in a different 
region, at the rear of this bow shock. On the other hand, [Oi], 
OH and high-yupC O line s correlate well with the [Feii] at 26 
fim (Neufeld et al ..l2009h. th e H2 ro-vibrational fine v =1-0 S(l) 
(Caratti o Garatti et al.L 12006!) tfiat traces fiot gas (T >2500 K) as 
well as witfi tfie pure H2 rotatio nal lines trac ing gas at tempera- 
tures between 200 and 1400 K dNisini et al.L |2010a). This indi- 
cates that the FIR PACS lines trace a more compact and more 
excited region tfian tfie bulk of tfie Bl bow sfiock, most prob- 
ably tfie "jet-sfiock" wfiere tfie jet powered by tfie protostar is 
impacting on tfie internal cavity wall. 



' HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, The NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 



3.2. Physical conditions In the hlgh-J^p CO emitting gas 

In order to look for possible spatial gradients of the physical 
conditions in the CO gas, we calculated the excitation tempera- 
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Fig. 2. Left: PACS map of CO (14-13) in LI 157-Bl (gray scale) 
compared with SiO(2-l) (red contours) and CH3CN (8^—7 a:) 
(green contours). Right: PACS map of [Oi]63/im(gray scale, 
first level is 3x10 W m^^/spaxel, corresponding to 3 cr, steps 
are 7x10"'^ W m^^spaxel) compared with W2V =1-0 S(l) (red 
contours). Crosses mark the centre of the spaxels of the PACS 
FOV. 



tures, which gives a lower limit to the kinetic temperature if the 
gas is not in local thermodynamic equilibrium (LTE). To correct 
for the different beam sizes, we convolved all CO line maps to 
the resolution of the transition with the longest wavelength, i.e. 
13'.'1 (PSF at 186 jim) and then extracted the fluxes at the posi- 
tion of each spaxel. The fluxes for spaxels (-5", 7") and (0",0") 
are reported in Table I A. ll We then constructed the correspond- 
ing rotation diagrams for each spaxel, as illustrated in Fig.[3]for 
three positions. The CO lines in Bl are optically thin for J >9 , 
as in dicated by the '^CO lines detected with HIFI dLefloch et al.L 
12012 ). thus the slope of the diagram gives directly 1/rex. We find 
a similar excitation temperature Tex - 210 ± 10 K at all posi- 
tions, indicating that the excitation conditions are fairly uniform, 
although at off'set (9",2") the excitation temperature is slightly 
higher, 241 K. In fact, a secondary peak of H7 v =1-0 S(l) i s 
observed at this position (see Fig.|2l where lNisini et al.l ( l2010ah . 
from the H2 lines, derived a temperature about 100 K higher 
than at the CO peak. The total column density of CO averaged 
over 13'.'1, assuming LTE, ranges from 2x10'^ cm"^ at (9",2") 
to 6x10^5 cm-2 at (-5",7"). 

Constraints on the kinetic temperature and density of the 
CO gas were obtained us ing the radiative transfer code RADEX 
(IVan der Tak et al.Ll200l in a plane-parallel geometry. We used 
the collisional rate coefficients with H2 of ' Yang et alJ (1201 Ol) 
and built a grid of models with density between 10"* cm"^ and 
10** cm"^ and temperature between 10 K and 10"* K. In Fig.|4] 
the region of minimum is shown as a function of density and 
temperature for the CO peak at (-5",7"). The best fit is obtained 
forn =2x10^* cm"-' and rkin=1080 K, but other fits at lower tem- 
peratures combined with higher densities are also possible (see 
Fig. lA.4b . including the LTE solution. 

Additional constraints on physical conditions are provided 
by the CO lines observed with HIFI towards the (0",0") position 
that, as shown by the rotational diagram, has the same excitation 
condition of offset (-5", 7"). Indeed, line profiles in CO(16- 
15) and CO (13-12) reveal a specific spectral sh ape for the ex- 
cited CO component dominating the PACS range dLefloch et all 
I2OI2I) . Its flux contribution to CO (10-9) can then be extracted 
and scaled to 13'.'1 by assuming the same spatial distribution 
as for CO (16-15). The simultaneous fitting of the HIFI and 
PACS fluxes clearly excludes low-density / high-temperature so- 
lutions (Fig.|4]right panel). The best fit now has 7']iin=580 K and 



offset (-5", 7") T„=214 K N(C0) = 6xi0'= err 
offset (0",0") T„ = 207 K N(C0) = 5x10'^ cm" 
offset (9", 2") T„=241 K N(C0) = 2x10'^ cm" 
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Fig. 3. CO rotational diagrams at 13" resolution, with superim- 
posed linear fits. Diamond, squares and triangles correspond to 
off'set (0",0"), (-5",7") and (9",2"), respectively. The infeiTed 
excitation temperatures and total CO column densities are indi- 
cated. 
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Fig. 4. X distribution at the CO peak (left panel) considering 
only PACS lines and at the FOV centre (right panel) consider- 
ing PACS and HIFI lines. Contours indicate 68, 90, and 99 % 
confidence levels. 



«=2xl0^ cm ^ but still with a wide range of acceptable param- 
eters, 200< Tkin <1 100 K for n > 6 x lO^cm"^ (Fig.|A31l- 

Tighter limits come from the comparison of the CO column 
density at (-5",7") inferred from RADEX fits with the H2 col- 
umn density at the same position and temperature derived from 
Spitzer data (Nisini et al., 2010a), averaged over the same beam. 
We find that the ratio A?(CO)/A?(H2) is (0.3-2)xl0"'* for the best 
CO fits at Tkin of 200 and 580 K, respectively. In contrast, the 
CO relative abundance increases to unphysically high values for 
Tkin > 800 K. This clearly favours dense and moderately warm 
gas (200< Tkin <800 K, and n > 10^ cm"^). From these best-fit 
models we calculate a total cooling of L(CO) = 0.004 L© at the 
CO peak and an intrinsic CO column density A^(CO)= 2x10^^ 
cm"^, assuming the size of 7" derived from PACS maps. 

3.3. [0\] and OH lines: Evidence for a dissociative jet-sliock 

In the left panel of Fig. |5] the surface brightness of the 
[Oi]145^(m line (obtained using an upper limit to the size 
of 7") is compared wit h sho ck-model predictions from 
iFIower & Pineau des Foreli (12010) . The [Oi] 145/63 Une ratio 
indicates that both lines are optically thin, but we used the 
[O i] 145yum line for this comparison, because the 63jL/m line may 
be affected by for eground absorption, as usually found by ISO 
dLiseau et all 12006 ). Evidently only dissociative J-type shocks, 
with velocity Vs > 30 km s"' for n=2xl0'* cm"^ and Vs > 20 
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Fig. 5. Comparison between observed value (horizontal dotted 
line) a nd theoretical values for the Flower & Pineau des ForetH 
(I2OIOI) shock models (connected symbols) of the [Oi]145//m 
surface brightness (left panel) and the cooling ratio 
L([Oi]145+63)/L(CO) (right panel), as a function of shock 
speed. Solid lines refer to J-shocks and dashed lines to C- 
shocks; squares refer to pre-shock density «=2xl0'* cm"^ and 
triangles to «=2xl0^ cm"-'. 



km s~' for n=2xl0^ cm"-', can match the observed [Oi]145yum 
brightness. Non-dissociative shocks (of both J- and C-type) un- 
derpredict this line by at least two orders of magnitude. This 
conclusion is also supported by the detection of [Fe 11] at 26 fim 
and [Si] at 25 yum (Neuf eld et aU ,2009) at the [Oi] peak; the 
iron and sulfur fine structure lines are excited only in dissocia- 
tive shocks. In particular, the [Oi]/ [Si] observed intensity ra- 
tio is consistent with both lines originatin g from a J-shock with 
the same conditions as inferred above (Hollenbach & McKee, 
11989.) . The same models underpredict the [Oi]/[Feii] but this 
is likely because the iron is partially depleted on grains. Shock 
models reproducing the [O i] brightness predict a ratio of CO to 
[O i] total luminosity within a factor of a few of the observed one 
(see Fig.|5]right panel), therefore the dissociative J-shock traced 
by [O i] may also contribute to the high-7up CO emission. 

The sum of the three OH doublets at 79 yum, 84 //m, and 
119 fim (see Table IA.2l i amounts to three times the [Oi]145yum 
line, c onsistent with both species tra cing the same dissociative J- 
shock jHoUenbach & McKej,ll989l) . Moreover, the ratio of OH 
to CO lumino sity is an order of m agnitude higher than predicted 
for C-shocks (iKaufman & Neuf eld. 1996). Indeed, in C-shock 
all oxygen not locked in CO is converted in water, while in dis- 
sociative J-shocks the presence of FUV photons produced at the 
shock front delays the full conversion of free oxygen into water 
and maintains a significant abundance of OH. Interestingly, both 
the relative intensities of the OH doublets in L1157-B1 and the 
ratios of each OH doublet to [Oi]145/im and to Lboi = SL©, are 
very similar to those observed much closer to lo w-mass proto- 
stars dWampfler et al l l201(]l:lHerczeg et al.Ll20Tll) . This suggests 
that, OH may not be excited by infrared pumping, but predom- 
inantly by dissociative J-shocks in outflows far away as well as 
close to the source. 

In dissociative J-type shocks, the [Oi]63jum luminosity is 
proportional to the mass flux into the shock, whic h is close to 
the w ind mass-flux in a slow moving wind-shock ( HollenbachL 
Il985h . Considering the total flux of the [Oi]63yum line in the 
PACS FOV (10"'^ W/m-2), we derive L([Oi]63//m) = 0.002 
Lq and Mw([Oi])= 2x10"^ Moyr"' in the blue lobe. We may 



compare this value to independent estimate derived from the 
momentum and age of the blue outflow lobe traced by CO(l- 
0) by Bachiller et al. (2001). Rescaling the value these authors 
derived for a distance of 440 pc to our adopted distance of 
250 pc, we infer M„(CO) - 6x10"^ Moyr"' under the assump- 
tion of outflow-wind momentum conservation and wind veloc- 
ity of 100 km s"'. Given the large uncertainties in the outflow 
age (Downes & Cabrit, 2007) and possible self-absorption in the 
[Oi]63//m line, the two estimates are consistent and agree rea- 
sonably well with a jet-driven shock accelerating the outflow as 
the main excitation mechanism of [O i] in B 1 . 

4. Conclusions 

In the Herschel-PACS spectrum of the Bl shocked region of the 
LI 157 outflow, emission lines from CO, H2O, OH and [O i] have 
been detected. The analysis of the high-7up CO lines with a radia- 
tive transfer code has allowed us to constrain the physical condi- 
tions of the warm (200< Tkin < 800 K) and dense (n > 10^ cm"^) 
gas. The spatial resolution of PACS, unprecedented at FIR wave- 
lengths, has shown that CO, OH and [O i] all surprisingly peak at 
the rear of the bow shock, where others tracers of J-type shocks 
such as ro-vibrational H2 and [Fe 11] were also detected. The in- 
tensity of the bright [O i] lines clearly indicates for the first time 
that a dissociative "reverse" wind-shock is the dominant exciting 
mechanism of these lines. The L(CO) and L(OH) are also con- 
sistent with a J-type shock. Our analysis reveals that the C-shock 
previously found in B 1 appears to give only a minor contribution 
to the emission of [OI], OH, and high-ZupCO lines detected in 
the PACS spectra. 
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Fig. A.l. Spectrum of the spatial pixel at the centre of the FOV, 
RA(J2000) = 20*'39'"10:2, Dec(J2000) = +68°01'10'.'5. 

Table A.l. List of CO line fluxes at off^set (0",0") and (-5",7"), 
derived from the line maps convolved at 13'.'1 resolution. 



Line 




Flux offset (-5", 7") 


Flux offset (0",0") 






(10-=' Wcm-2) 


(10^=' Wcm-2) 


CO (22-21) 


118.58 


3.2+0.5 




CO (21-20) 


124.20 


4.2+0.9 


2.9+0.7 


CO (20-19) 


130.39 


4.3+0.6 


3.0+0.5 


CO (19-18) 


137.23 


5.4+0.7 


3.9+0.8 


CO (18-17) 


144.79 


6.3+0.8 


4.5+0.6 


CO (17-16) 


153.28 


9.0+1.0 


5.3+0.7 


CO (16-15) 


162.82 


12.2+1.4 


7.6+1.0 


CO (15-14) 


173.63 


13.4+1.6 


10.0+1.2 


CO (14-13) 


186.01 


15.9+1.7 


10.3+1.1 


CO (13-12) 


220.25 




10.0+1.6 


CO (10-9) 


260.24 




6.1+1.0 



Appendix A: 



In Fig. lA.ll and IA.2I we show the PACS spectra at the cen- 
ti-e of the FOV (the nominal position of Bl) RA(J2000) = 
20^39™10^;2, Dec(J2000) = +68°01'l(y.'5 and towards the PACS 
emission peak, at position RA(J2000) = 20*'39"9!3, Dec(J2000) 
= -)-68°01'17'.'8, located at off'set (-5"; 7") with respect to the 
centre of the FOV. 

In TableOwe list the fluxes at off^set (0",0") and (-5",7") 
extracted from the line maps convolved at the diffraction size 
of the CO transition at the highest wavelength, i.e. 13'.'1 (PSF 
at 186 fim). They were obtained by multiplying the line surface 
brightness by an area of 1.133 x 13'.'!^, and were used for the ro- 
tational diagram and the RADEX modelling, assuming the same 
source area. 
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Fig. A.2. Spectrum of the spatial pixel centred at RA(J2000) = 
20'^39'"9!3, Dec(J2000) = -h68°01'17'.'8, coiTesponding to an off- 
set (-5",7"). 

Table A.l. OH hne fluxes measured at offset (-5",7") after PSF 
correction for a point source. The FWHM was set to the instru- 
mental value. 



Line 




FWHM 

(yum) 


Flux 
(10-"W/m=) 


OH^Hi 




1 - 
2 


3+ 
2 


79.11 


0.039 


3.5+0.9 


OH^ni 




1 + 

2 


3- 
2 


79.18 


0.039 


3.1+0.9 


OH ^n, 


-^n^y = 


7+ 
2 


_ 5- 


84.42 


0.038 


3.1+0.9 


OH ^Hj 

-> 


--nij = 


7- 
2 


3+ 
2 


84.58 


0.038 


3.4+0.9 


OH^n^ 


-^n\j = 


5- 
2 


3 + 
2 


119.23 


0.117 


3.0+0.4 


OH^Hj 

-> 




5+ 
2 


3- 
2 


119.44 


0.117 


4.2+0.4 


OH^ni 




3+ 
2 


1 - 
' 2 


163.12 


0.126 


< 1.2 


OH^Hi 


-^Hiy = 

1 


3- 
2 


1 + 
2 


163.40 


0.126 


< 1.2 



In Table IAT2I we list the fluxes of the OH lines measured at 
offset (-5",7"). Since the OH lines were detected only in one 
spaxel, we applied the flux correction for a point source that 
takes into account the size of the PSF. 

In Fig. IA.3I we show the comparison between the CO line 
fluxes at the centre of the PACS FOV, off'set (0",0"), and three 
models: the best-fit model obtained considering only the PACS 
lines (/up >14) with n =2x10"* cm"^ and T^^ =1080 K; the best- 
fit model obtained considering also the HIFI CO (13-12) and 
(10-9) fines with n =2x10^ cm"^ and 7^,, =580 K and the LTE 
model with kinetic temperature equal to the excitation tempera- 
ture derived from the rotational diagram, with n =2x10^ cm"-' 
and Tkin =207 K. In Fig. lA.4l the same models are shown com- 
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Fig.A.3. Fit of the CO lines at offset (0",0"), taking into ac- 
count only PACS data (dashed curve), adding the HIFI lines 
(solid curve) and assuming the LTE temperature (dotted curve). 
Squares are the PACS fluxes and bars represent the range of 
fluxes estimated from HIFI lines. Error bars include the statis- 
tical error from the line fitting and the 10 % calibration uncer- 
tainty. 




Fig.A.4.Fit ofthe CO lines at off'set (-5",7"). Error bars include 
the statistical error from the line fitting and the 10 % calibration 
uncertainty. 



pared to the CO line fluxes at the peak of the CO emission, i.e. 
ofi^set (-5",7"). 



